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ABSTRACT 

It is shown that the radio and gamma-ray emission observed from newly-found "GeV-bright" supernova 
remnants (SNRs) can be explained by a model, in which a shocked cloud and shock-accelerated cosmic rays 
(CRs) frozen in it are simultaneously compressed by the supernova blastwave as a result of formation of a 
radiative cloud shock. Simple reacceleration of pre-existing CRs is generally sufficient to power the observed 
gamma-ray emission through the decays of 7r°-mesons produced in hadronic interactions between high-energy 
protons (nuclei) and gas in the compressed-cloud layer. This model provides a natural account of the observed 
synchrotron radiation in SNRs W5 1C, W44 and IC 443 with flat radio spectral index, which can be ascribed to 
a combination of secondary and reaccelerated electrons and positrons. 

Subject headings: cosmic rays — acceleration of particles — ISM: individual objects (W51C,W44,IC 443) — 
radiation mechanisms: non-thermal 



1. INTRODUCTION 

Luminous extended GeV 7-ray emission associated with 
middle-aged supernova remnants (SNRs) has recently been 
unveiled by the Large Area Telescope (LAT) onboard the 
Fermi Gamma-ray Space Telescope. Specifically, SNRs 
W51C, W4 4, IC 443, and W28 are sp atially resolved with the 
Fermi LAT (lAbdo et al.ll2009ll2010allblkl) . The four SNRs are 
interacting with molecular clouds, as evidenced e.g., by 1720 
MHz OH maser emission. The radio and 7-ray emission from 
these SNRs share similar characteristics. The synchrotron ra- 
dio emission has a large flux of 160-310 Jy at 1 GHz with flat 
spectral index of a ~ 0.26-0.40. The GeV 7-ray spectrum 
commonly exhibits a spectral break at around 1-10 GeV, and 
the luminosity ranges L 7 = (0.8-9) x 10 35 erg s" 1 in the 1-100 
GeV band. Other cloud-interating SNRs associated with the 
LAT sources, such as CTB 37A, also emit 7-ray s at a lumi- 
nosity of L 7 - 10 35 erg s" 1 (ICastro & Slandl2010l) . 

A prototypical example of the GeV-bright SNRs is SNR 
W44. The radio continuum map of W44 exhibits filamen- 
tary and sheet-like structures of synchrotr on radiation well 
correlated with the s hocked H2 emissio n dReach et al.| l2005t 
ICastelletti et al.1 12007b . According to iReach et al.l (120051) . 
the bulk of the synchrotron radiation can be ascribed to a 
fast molecular shock of a velocity v s ~ 100 km s" 1 advanc- 
ing through a molecular cloud of a preshock density no ~ 
200 cm -3 . By passage of the blastwave of the supernova rem- 
nant, the shocked molecular cloud forms a thin sheet due to 
radiative cooling. The radio filaments are thought to come 
from the compressed zone behind the shock front. 

The synchrotron radio emiss ion arising from such a 
"crushed cloud" was modeled by Bland ford & Cowiel (1982). 
It was shown that reacceleration of pre-existing CR electrons 
at a cloud shock and subsequent adiabatic compression results 
in enhanced synchrotron radiation, c apable of explaini ng the 
radio intensity from evolved SNRs. iBvkov et al.l d2000l) dis- 
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cussed a similar scenario, in which direct electron accelera- 
tion from the thermal pool is also invoked. 

The 7r°-decay 7-ray emission should be enhanced in the 
crushed clouds in the same manner as the synchrotron radi- 
ation. In this Letter, we demonstrate that the newly-found 
7-ray emission from middle-aged SNRs can be readily un- 
derstood within the crushed cloud scenario, in which the 7- 
ray emission comes from shocked clouds overrun by SNR 
blastwaves. It should be emphasized that our model is essen- 
tially different from the scenarios adopted i n the recent papers 
(e.g., iFuiita et all 120 lot [Torres et al] |2010l) . where molecular 
clouds in the vicinit y of SNRs are assumed to be illuminated 
by ru naway CRs (lAharonian & Atoyanl Il996t iGabici et"aTJ 
2009). While the GeV and TeV 7-rays outside the southern 
boun dary of SNR W28 (lAharonian et al] |2008t iGiuliani et ail 



2010; lAbdo et al.ll20l0c¥ may be explained by such runaway 
CRs, we argue here that the luminous GeV 7-ray emission 
in the directions of cloud-interacting SNRs emerges from the 
radiatively-compressed clouds. 

2. THE MODEL 

2.1. Cloud Shock Structure 

Let us consider a strong shock driven into a molecular cloud 
by the high pressure behind a supernova blastwave. Using the 
number density of hydrogen nucleus in the preshock cloud, 
no, the preshock magnetic field is described by a dimension- 
less parameter b: 



Bq = b\ (no/cm 3 ) /zG. 



(1) 



Zeeman measurements of self-gravitating molecular clouds 
show that b is roughly cons tant from one cloud to another 
with b ~ 1 (Crutcheri ll999l) . We are concerned primarily 
with a fast (v s > 50 km s" 1 ) /-type shock (Drainel [T980l) . in 
which ambipolar diffusion and radiative cooling are unim- 
portant in a shock dissipation layer. The steady-state post- 
shock structure of a fast m olecular shock is described in 
iHollenbach & McKed dl989l) . which is referred to as HM89 
hereafter. 

The initial temperature immediately behind the shock front 
is T s ~ 3.2 x 10 s v^ 7 /x,K, where v S 7 = v s /(100km s" 1 ), andx, 
is the number of particles per hydrogen nucleus (HM89). The 
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density rises from a preshock density «o by a factor of r s h = 4, 
in the strong shock limit. For v S 7 > 1.2, ionizing radiation 
produced in the immediate postshock layer is strong enough 
to fully predissociate and preionize the upstream cloud. 

As the gas radiatively cools downstream, the temperature 
decreases and the density increases. The compression of the 
cooling gas is limited by the magnetic pressure in the cases 
that we shall be considering. The density of the cooled gas, 
n m , is determined as 

n m ~ 94«oVs7^ _I , (2) 

which is obtained by equating B 2 /87r with the shock ram 

pressure nopnv 2 , where B„, = ■y/2/3(n,„/«o)-Bo is the com- 
pressed magnetic field, and /jh is the mass per hydrogen nu- 
cleus (HM89). (We have assumed that the preshock magnetic 
field is randomly directed.) The [O I] line (63/zm) surface 
brightness of a face-on /-shock is proportional to the parti- 
cle flux into the shock, hqv s . The peak surface brightness 
of the [O II line is 1 x 10" 3 erg cm" 2 s" 1 sr" 1 in SNR W44 
([Reach &Rholll 996). and a factor of 2 smaller in IC 443 
dRho et al.ll2001l) . indicating the particle flux of order nov s ~ 
10 9 crrrV 1 (HM89) and therefore n m ~ 9 X 10 3 /?" 1 cm" 3 . 

A certain column density, 7V coo i, has to be transmitted by a 
shock for gas to cool down to 10 4 K and become radiative. 
This occurs in a column density N c „„\ ~ 3 x 10 17 v 4 7 cm -2 for 
v,7 = 0.6-1.5 dMcKee et al.lll987l) . Recombination and pho- 
toionization by the ultraviolet radiation produced upstream 
are balanced in 7V c00 i < N < Mon ~ 10 19 cm -2 . Beyond Ni 0D , 
the ionizing photons are absorbed and molecular chemistry 
commences. By setting a typical elapsed time since shocked 
to be t c = t/2, where t = 10 4 ?4yr is the age of the rem- 
nant, and using «o.2 = "o/(100 cm" 2 ), the column density 
of the compressed cloud is written as N c = nov s t c ~ 1.5 x 
10 20 no,2V.,7f4 cm" 2 . 

The bulk of synchrotron radio waves and 7-rays should be 
emitted in the compressed gas with a constant density n m 
(within a factor of 2) and magnetic field B m . Note that we 
idealize the shock as one-dimentional and ignore any effects, 
such as lateral compression, caused by the secondary shocks. 
The secondary shocks could affect the leakage of high-energy 
particles from the compressed cloud, playing an indirect role 
in the 7-ray production. Also, the ultraviolet radiation pro- 
duced in the secondary shocks may change the ionization 
level of the precursor of the main cloud shock. 

2.2. Shock Acceleration and Adiabatic Compression 

Following Blandford & Cowie ( 1982), we consider a con- 
servative case in which only pre-existing CRs are acceler- 
ated by the process of diffusive shock acceleration at a cloud 
shock. Suprathermal particles may be injected to the accel- 
eration process at the shock front, despite slow acceleration 
and fast Coulomb losses. However, we shall show below that 
the shock acceleration of pre-existing CRs alone appears to 
suffice as the origin of the observed 7-ray emission from the 
cloud-interacting SNRs, and therefore consider the reacceler- 
ation case only. 

Let n acc (p)dp be the CR number density in p ~ p + dp, 
transmitted by a shock, and «gcr(p) be the pre-existing am- 
bient CR de nsity. According to the theo ry of diffusive shock 
acceleration (Blandford & Eichler 1987), for p < Pbr/max (see 
below), 

«acc(p) = (a + 2) p~ a f dp'nacRip') p' (a ~'\ (3) 
Jo 



where a = (r s h + 2)/(r s h- 1) and r s h is the shock compression 
ratio, which is assumed to be r s h = 4. 

Assuming that the density of the Galactic CRs in the molec- 
ular cloud is same as that in the general interstellar medium, 
we adopt the Galactic CR proton spectrum of the form: 

«GCR. P (p) = 47r/ ;) /3 L5 /?o 2 ' 7f \ (4) 

where p = p/(GeV/c), J p = 1.9cm" 2 s" ! sr" 1 GeV" 1 , /3 is the 
proton velocity in units of c. A low-energy cutoff at a kinetic 
energy of 50 MeV is ap plied. This sp ectrum lies in between 
Shika ze et all (120071) and lStrong etaLl (12004 at 100 MeV. For 
the CR electron+positron spectrum, we use: 

nGCR, e (p) = 47r7 f ,po 2 (l+po) °' 55 , (5) 

extending down to a low-energy cutoff of 20 MeV at 
which point ionization losses in the Galaxy should make 
the spectrum fiat. The normalization factor is J e = 2 x 
lO^cm-V^r'GeV" 1 . 

The preshock upstre am gas is on l y part ially ionized for 
Vj7 < 1.2. Recently, Malkovetal. (2010) have proposed 
that strong ion-neutral collisions accompanying Alfven wave 
evanescence lead to steepening of the spectrum of accelerated 
particles; the slope of the particle momentum distribution be- 
comes steeper by one power above p\, T = 2eBoVA./cz4-n, where 
Va is the Alfven velocity and ^;_ n ~ 9 x 1Q~ 9 n n oT^ >A s" 1 is 
the ion-neutral collision frequency. Here n n .0 denotes the den- 
sity of neutrals in units of cm" 3 and T4 is the precursor tem- 
perature in units o f 10 4 K. Interestingly, the parameters of 
iReach et"afl (120051) estimated for the radio filaments of W44, 
namely v s i ~ 1 and «o.2 ~ 2, together with Bo = 30 pG and 
T4 = 1 predict pt,r ~ 10 GeV/c using the precursor ioniza- 
tion fraction calculated by HM89. This agrees with t he break 
value measured by the Fermi LAT (Abd o et al.l2010al) . We in- 
troduce spectral steepening to H acc (f>), by multiplying a factor 
°f Pbr/ P above p bl . 

Even if the preshock gas is fully ionized, spectral steep- 
ening due to a finite acceleration time would be unavoid- 
able already in the Fermi bandpass. The timescale of diffu- 
sive shock acceleration can be written as f acc ~ (10/3)ricr g v~ 2 , 
where r g = cp/eBo is the gyroradius and 77 > 1 is the gy- 
rofactor; r\ ~ 1 has been obtained in young SNRs like 
RX J1713.7-3946 (lUchiyamaet al .1120071) . Equating f acc with 
t c , the maximum attainable energy is obtained as cp m . dx ~ 
50(7 ? /10)" 1 v 2 v B- 5 f4 GeV, where B_ 5 = B /(10" 5 G). A fac- 
tor of exp[-(/?//? max )] is multiplied with n acc {p) to introduce 
the maximum attainable energy. 

The high-energy particles accelerated at the shock expe- 
rience further heating due to adiabatic compression, as the 
gas density increases until the pressure is magnetically sup- 
ported. Each particle gains energy as p — > s^ 3 p, where 
s = (« m /»n)/r s h, and the dens ity increases by a factor of s 
(see iBlandford & CowidlT982l) . Therefore, the number den- 
sity of accelerated and compressed CRs at the point where the 
density becomes ~ n m is 

«ad(p) = S 2/3 « acc (s- 1/ V). (6) 

2.3. Evolution in the Compressed Region 

In the compressed region of the cloud with a constant den- 
sity ~ n m , high-energy particles suffer from energy losses 
such as Coulomb/ionization losses. Also, the production of 
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TABLE 1 
Model Parameters for SNR W44 



Parameters 



Values 



Assumed SNR Dynamics 



Distance: D 
Radius: R 

Age: ( 

Explosion energy: B51 



2.9 kpc a 
12.5 pc a 
(6 = 15') 
10000 yr a 
5 a 



Preshock Cloud Parameters 

Density: «o 
Filling factor: / 
Magnetic field: Bo 



200 cm 41 
0.18 
25 fiG 



Dependent Parameters 



Cloud shock velocity: v s 
Break momentum (T4 = 2): pbr 
Maximum momentum (rj = 10): p m 



100kms-' a 
7GeV/c 
122 GeV/c 



a Taken from lReach et al.l d2005r) . 

secondaries in inelastic proton-proton collisions can be a sig- 
nificant source of high-energy electrons and positrons, since 
the energy loss timescale due to the pp collisions, t pp ~ 
6 x 10 7 (n m /cm" 3 ) -1 yr, generally becomes comparable to t c . 

Let N p (p,t)dp represent the number of protons in p ~ 
p + dp integrated over the emission volume at time t. We em- 
ploy a usual kinetic equation to obtain the proton spectrum 
N p (p,t): 



dN p (p,t) _ d 
dp 



dt 



[b(p)N p (p,t)]+Q p (p), 



(7) 



where Q p is the proton injection rate, and b(p) = —p represents 
the proton energy losses. Charged particles with p < p mdx 
are assumed to be effectively trapped within the compressed 
gas. The particle spectrum for accelerated electrons, N e (p,t), 
and for secondary e ± resulting from hadronic interactions, 
N sec (p,t), are computed similarly. 

The injection rates of primaries are related to n a d(p) : 



Q(p)= «ad(/>), 



Tlmfr 



(8) 



where V = 4irR 3 /3 is the SNR volume with a radius R, and 
/ is the preshock filling factor of the compressed cloud. The 
volume of the preshock cloud is then fV. For simplicity, we 
assume that Q(p) for primaries is time-independent; v s is as- 
sumed to be constant and p m . dx is evaluated simply at t = t c . 
The injection rate of secondary e^, Qsec(Pe>0> i s determined 
by N„( p D , t) following the prescription given by Kam ae et al.l 
(l2006h . 

2.4. SNR Dynamics 

It may be desirable to relate the cloud shock velocity with 
the physical parameters that describe the blastwave of the 
remnant. We consider a remnant in the Sedov staged 



^12.5 = 



£51 

«a.O 



1/5 



(9) 



4 This implies t <t lr , where f lr is an age for transition to the radiative phase 
of SNR evolution iBlo ndin et alTO^ : f tt ~ 3 X 10 4 E*{" fT 9 /' 7 yr. 



where R12.5 = R/(12.5 pc), £51 denotes the kinetic energy 
released by the supernova in units of 10 51 erg, and n a .o = 
« a /(cm" 3 ) represents the ambient (intercloud) density. The 
blastwave velocity is Vb = OAR/t. When a molecular cloud 
is struck by the blastwave, a strong shock is driven into the 
cloud with a shock velocity of 



v s = k 



no 



1/2 



;65 ^51 -^12.5 s 



(10) 



where k ~ 1.3 is adopted. N ote that k depends weakly on 
v s /vb (Mc Kee & Cowielll975l) . ranging k = 1—1.5 in the cir- 
cumstances of our interest. To drive a fast molecular shock of 
v S 7 > 0.5 into a molecular cloud, a condition of 



n 7<k 2 E 5l R 



12.5 



(11) 



must be satisfied. 




10 u 10 1 
E [GeV] 

Fig. 1. — 7r°-decay 7-ray spectra calculated for the reacceleration 
model using various sets of parameters. Spe ctral data points are 
for SNR W51C (fi lled circles: lAbdo et al.l2009l). W44 (open circles: 
lAbdo et alj20 10a). and IC 443 (squares: I Abdo et al. 2010b). 
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Fig. 2. — Radio (left) and 7-ray (right) spectra of SNR W44 together 
with the reacceleration mod el using the param eters in Table [T] most 
of which are copied from Reach et al. (2005). The radio fluxes are 
scaled by a factor of 0.5 (see the text). 
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3. RESULTS 

Let us calculate the nonthermal radiation arising from the 
radiatively compressed clouds. The gas density and mag- 
netic field strength are constant over the emission volume in 
which high-energy particles are distributed with the volume- 
integrated spectrum of N(p,t c ). Leptonic components include 
synchrotron radiation, inverse-Compton scattering, and rela- 
tivistic bremsstrahlung; the inverse-Compton emission is neg- 
ligible as compared to the bremsstrahlung component (see, 
e.g. jAbdo et al.ll2010al) . Inelastic collisions between CR pro- 
tons/nuclei and gas nuclei lead to 7r°-decay 7-ray emission, 
which constitutes the main component in the Fermi-LAT band 
under the assumption that only pre-existing CRs can be accel- 
erated at the cloud shock. 

We do not consider the densest part of the interacting 
molecular cloud. For example, det ections of OH(1720 MHz) 
maser emission dFrail et al.l [T996) indicate the presence of 
slow C-type shocks in dense molecular clumps, say v s j ^0.3 
and «o.2 ~ 10, in addition to the fast J- type shock described 
in ^2] However, the 7-ray emission from the slow non- 
dissociative shocks is not expected to be strong because of in- 
effic ient shock-accelera tion and weaker gas compression (see 
also iBykov et al.ll2.Q00h . Indeed, the contribution from such 
dense molecular cl umps to the total r adio intensity appears 
small in SNR W44 (iReach et al .1120051) . The largest preshock 
density we should consider is given by Eq. (fTTT i. Also, we ig- 
nore the blastwave region, since the radiative compression is 
essential in our model. 

3.1. Gamma-ray Luminosity 

We show here that the 7-ray luminosity anticipated within 
this scenario agrees well with the observed luminosity of 
~ 10 35 erg s" 1 . In Figure Q] the 7-ray spectra of SNRs W5 1C, 
W44, and IC 443 measured with the Fermi-LAT are shown 
in the so-called vL u form in un its of 10 3S e r g s" 1 . The LA T 
spectral points are taken from lAbdo et al.l d2009l l2010allbh 
and converted into the vh v form using the distances of 6 kpc 
(W51C), 2.9 kpc (W44), and 1.5 kpc (IC 443). 

To demonstrate the expected level of the 7-ray luminos- 
ity, we present the spectra of 7r°-decay 7-rays with varying 
R, no, and £51. The following parameters are fixed: b = 2, 
/ = 0.2, n a fi = 1. Also, p max is set by adopting r\ = 10, and 
Pbr is set by T4 = 2 and by the ionization fraction calculated 
based on HM89, here and hereafter. The 7-ray spectra sim- 
ply scale as oc /, and depend very weakly on n a . The black 
lines in Fig. Q] show the results obtained for »o,2 = 0.3 (solid 
curve) and «o.2 = 3 (dashed curve) in the case of R = 10 pc 
and £51 = 1. The 7-ray luminosity around 1 GeV varies only 
within a factor of ~ 2 between «o.2 = 0.3 and «o.2 = 3, while 
that at 100 GeV changes more than an order of magnitude. 
On the other hand, the blue lines in Fig. [T] show the spectra 
calculated for R = 5 pc (solid curve) and R = 15 pc (dashed 
curve) in the case of «o,2 = 1 and £51 = 1. The 7-ray luminos- 
ity differs by a factor of ~ 3. We note that the shocked cloud 
mass amounts to ~ 10 4 M© in the case of R = 15 pc. Finally, 
to explore the most luminous scenario, we adopt £51 = 5 to- 
gether with R = 30 pc and «o,2 = 1 (red curve). The 7-ray 
luminosity reaches ~ 10 36 erg s" 1 , in good agreement with 
the observations of SNR W51C, which is indeed the most lu- 
minous SNRs in gamma-rays. Our model generally predicts 

£ 7 ~ 10 35 {f /Q.2)El[ 3 erg s" 1 to the first order, which led us to 
conclude that the £en«/-detected 7-rays are quite likely due to 



the decays of 7r -mesons produced by the pre-existing CRs ac- 
celerated and subsequently compressed in the shocked cloud. 

3.2. Flat Radio Spectra 

The GeV-bright SNRs W51C, W44, and IC 443 are also 
radio-bright objects. As shell-type SNRs, their radio spectra 
are remarkably flat with spectral inde x of a ~ 0.26 (W51C : 
iMoon & Koo1ll994l) . a ^ 0.37 (W44: ICastelletti et al.l l2007l). 
and a ~ 0.36 (IC 443: lErickson & Mahonevl 1 19851) . with a 
typical uncertainty of 0.02, being inconsistent with a = 0.5 
that is expected by shock-acceleration theory. Our model nat- 
urally explains the flat radio spectrum. Let us demonstrate by 
presenting the radio and 7-ray modeling of SNR W44 how the 
radio an d 7-ray spectra can be simultaneo usly reproduced. 

Radio (ICastelletti et al.[2 007) and 7-ray (lAbdo et al.l2010al) 
spectra of SNR W44 are shown in Figure [2] One half of the 
total synchrotron flux measured for W44 is assumed to orig- 
inate in a fast molecular shock, wh ich is roughly consistent 
with Table 2 of Reach et al. (2005) . The rest is at tributed to 
the blastwave region. Table 2 of iReach et all £2005) was cho- 
sen as an initial set of model parameters: R12.5 = 1, h = 1, 
£gi=5, and «o,2 = 2. We then attempted to reproduce the non- 
thermal radiation spectra by varying / and Bq, and found that 
/ = 0.18 and Bo = 25 /iG provide a good fit to the data (see 
Table [T|and Fig. |2). The radio measurements can be recon- 
ciled with this model, in which the synchrotron radiation is 
largely contributed by secondary electrons and positrons. The 
flat radio spectrum is generically expected in our model. 

4. DISCUSSION 

The radiatively-compressed cloud provides a simple expla- 
nation for the radio and 7-ray data. Interestingly, the observed 
steepening in the 7-ray spectra is successfully reproduced by 
/?br in the case of SNR W44. However, there may be other 
explanations for the steepening. For example, high-energy 
particles may be prone to escape from the compressed mag- 
netized cloud. Also, the spectral break may be due to the fact 
that the crushed clouds have a range of no- A superposition 
of 7-ray spectra characterized by different /j max (as a result 
of different «o) could look like a break. We are primarily in- 
terested in understanding the 7-ray luminosity rather than the 
spectral shape in this paper, and therefore we did not explore 
this issue. 

The simple re-acceleration of pre-existing CRs and sub- 
sequent compression alone would not fully explain the 7- 
rays associated with cloud-interacting SNRs. For exam- 
ple, the GeV-TeV 7-ra y emission found outside the radio 
boundary of SNR W28 (lAharonian et al.ll2008b iGiuliani et all 
l2010t [Abdo et an i2010cl) may r epresent the molecular cloud 
illuminated by ru naway CRs (lAharonian & Atovanl 119961: 
iGabici et al .1120091) . Also, we assumed pre-existing CRs in 
the cloud to have the same spectra as the galactic CRs in the 
vicinity of the solar system. However, the ambient CRs in the 
pre-shock cloud may deviate from the galactic pool due to the 
runaway CRs that have escaped from SNR shocks at earlier 
epochs. If this is the case, modeling of the 7-ray spectrum, at 
TeV energies in particular, should take into account modified 
pre-existing CRs in the pre-shock cloud. 

We acknowledge the useful suggestions of the anonymous 
referee, which improved the manuscript. We wish to thank 
Heinz Volk and Felix Aharonian for valuable discussions. 
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